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Abstract 

The quark delocalization and color screening model, a quark potential 
model, is used for a systematic search of dibaryon candidates in the u, d and s 
three flavor world. Color screening which appears in unquenched lattice gauge 
calculations and quark delocalization (which is similar to electron delocaliza- 
tion in molecular physics) are both included. Flavor symmetry breaking and 
channel coupling effects are studied. The model is constrained not only by 
baryon ground state properties but also by the A^-A^ scattering phase shifts. 
The deuteron and zero energy di-nucleon resonance are both reproduced quali- 
tatively. The model predicts two extreme types of dibaryonic systems: "molec- 
ular" like the deuteron, and highly delocalized six-quark systems among which 
only a few narrow dibaryon resonances occur in the u, d and s three flavor 
world. Possible high spin dibaryon resonances are emphasized. 



I. INTRODUCTION 



Quantum chromodynamics(QCD) is believed to be the fundamental theory of 
the strong interaction. High energy processes are calculable due to the asymptotic 
freedom property of QCD. The majority of low energy processes are uncalculable due 
to infrared confinement. Lattice gauge calculations may suffice in the confinement 
regime, but will still suffer from large numerical uncertainties for the prediction of 
many hadron properties in the near future. This leads to a reliance on QCD inspired 
models to explore hadron physics for the time being and perhaps even in the future, 
due to the complexity of QCD. The existing models (potential, bag, soliton, etc.) are 
quite successful in understanding hadron (meson and baryon) properties, but have not 
been very successful for hadron interactions. Only recently have there been positive 
indications for obtaining the whole N-N interaction from QCD models'^'^l. 

An outstanding problem is the fact that all of these models, including lattice gauge 
calculations, predict that there should be: multiquark systems {qqy, q'^q, q^; quark 
gluon hybrids qqg, q^g; and glueballs in addition to the qq mesons and q^ baryons. 
Experimentally there are no well established candidates for these exotics. In a rela- 
tivistic theory, since quark and gluon number is not conserved, any meson state can 
be a mixture of qq, {qqy, and qqg; any baryon state can be a mixture of q^, q'^q and 
q^g. It is quite possible that these exotics, {qq^, qqg, g"^, q^q, q^g, exist in the normal 
meson and baryon states'^'. Polarized lepton nucleon scattering measurements have 
aroused a new round of hadron structure studies, wherein these exotic components 
are explored in connection with the normal q^ and qq components'^^ . However, q^ is 
really a new quark system sector, different from that of mesons and baryons. We call 
a baryon number B=2 state, which is quasi-stable, a dibaryon. Its minimum configu- 
ration is q^. Since Jaffe predicted the first dibaryon, the H particle'^l, a large number 
of dibaryon calculations with all the above-mentioned QCD models have been carried 
out and almost all support the existence of dibaryons'^' . If the present absence of an 
experimental dibaryon signal continues, then all these QCD models (and even QCD 
itself) should be questioned. Therefore, the dibaryon is a good place to test QCD 
and its models. 

Recently, Silvestre-Brac et al. reported a new systematic dibaryon calculation 
based on the chromomagnetic model'^l As pointed out by Lichtenberg and Roncaglia, 
the chromomagnetic model Hamiltonian is oversimplified t^' . The chromomagnetic in- 
teraction can give only the N-N short range repulsion but not any N-N attraction. 
Many dibaryon model calculations have the same deficiencies unless a phenomenolog- 
ical meson exchange is invoked. To study the dibaryon, it is better to have a model 
Hamiltonian which can at least fit the N-N interaction qualitatively. Then we can 
expect that such a model prediction may be relevant to real dibaryon states. Another 



1 



deficiency of many prevailing model calculations is that the model Hilbert space is 
rather restricted. In some model calculations^^'^'^^, the six quarks are assumed to 
be completely merged into a single confinement region (which we term a 'fully de- 
confined' model). In other model calculations, the quarks are assumed to be always 
confined separately in two distinct baryons (which we term a fully confined model) '^1. 
The real situation is quite possibly neither completely merged nor always separately 
confined, but rather, in between, i.e., partially deconfined due to the interaction of 
two baryons. A more realistic model calculation allows the six quark system to choose 
the preferred configuration by its own dynamics. 

To remedy these model deficiencies, we developed a model which we will call the 
quark delocalization, color screening model(QDCSM)P^. The model Hilbert space is 
enlarged to include the fully confined and fully deconfined models as two extremes and 
the real configuration is determined variationally by the dynamics of the six quark 
system. In this way the system is allowed to develop its own preferred distortion. 
The model Hamiltonian is sufficiently realistic to produce a qualitatively correct A^- 
N phase shift. We especially take into account the possible difference of the q-q 
interaction inside a hadron and between two colorless hadrons due to the nonlinearity 
of QCD (see section II). We use this model to do a systematic search within the 
u, d, s quark three flavor world, expecting it to improve the reliability of estimates 
on promising dibaryon candidates. This expectation is realized in the A^-A^ channels 
(see section IV). Of course we cannot expect that the model estimate is quantitatively 
correct, because both the model Hamiltonian and Hilbert space are restricted to be 
simple enough to do a systematic search. As emphasized by Silvestre-Brac'^^l this 
kind of systematic search serves the purpose of delimiting, among the thousands of 
multiquark states, the most promising candidates. Our intent is to assist experimental 
efforts to explore a challenging question in hadron physics by providing more reliable 
theoretical estimates. 

This paper is organized as follows. In section II, the model Hamiltonian and 
Hilbert space are described. Section III is devoted to a sketch of the calculation 
method. (A more complete description will be reported separately). The results are 
given in section IV and a conclusion in section V. 

II. QUARK DELOCALIZATION, 
COLOR SCREENING MODEL 

Quark potential models are quite successful in describing single hadrons; therefore 
we adopt the usual potential model Hamiltonian to describe a single baryon: 

3 2 3 

H{3) = ^(m, + ^) + E ViJ - Tc, (1) 

i=l ^"^i i<j=l 
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V ■ — V-^ + 

V"j = -Xi-Xjar^j, 

■' 4 \rij 2 \mf rrij Smiirij j j ^ y 

The symbols in eq.(l) have their usual meaning. For the confinement potential V^, 
we assume a quadratic form to simplify the calculation. A possible constant part is 
omitted to reduce the number of parameters. In the effective one gluon exchange 
potential V'^ ^ only the color Coulomb and color magnetic terms are retained, because 
we are only interested in the ground state. The effect of the momentum dependent 
Darwin term has been checked and found not to be critical; hence it is also omitted 
to reduce the calculational burden. 



A variational three quark wave function (WF) of the form 

z^(123) = x(123)r/5/j(123)0(123), (2) 

is assumed to describe the ground state baryons. Here x(123) is the color singlet WF, 
?75/j(123) is the symmetric spin-flavor SU2X/ DSU/XSU2 WF [S = strangeness, / = 
isospin, J = spin). 

0(123) = (l>is{n)Mr2)Mrs): (3a) 



Mr1=(^j e-^, (3b) 

s is a reference center, and 6 is a baryon size parameter to be determined by the 
stability condition 

aM(123) 



db 



= 0, (4) 



where M(123) = (V' ( 123) |//( 123) 1^^(123)). 

The other model parameters arc fixed as follows: the u, d quark mass difference 
is neglected and m = m„ = is assumed to be exactly ^ of the nucleon mass M, 
i.e., m — 313 MeV. The quark gluon coupling constant as is determined by the N- 
A mass difference. The confinement potential strength a is determined by the zero 
nucleon binding. The strange quark mass is determined by an overall fit to the 
strange baryon masses under the assumption that all the fiavor octet and decuplet 
baryons have the same r.m.s radius h. (Choosing different values of h for different 
baryons will make the calculation more elaborate and is left as a future refinement.) 
The fitted parameters are: m — 313 MeV, = 634 MeV, h — 0.603 fm, q;s=1.54, 
a = 25.13 MeV-fm~^. The theoretical baryon masses are compared with experimental 
values in Table I. 
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Table I. Baryon masses (MeV) 





N 


A 


A 


E 




S* 




Q 


exp. 


939 


1232 


1115 


1193 


1318 


1385 


1533 


1672 


theor. 


939 


1232 


1118 


1217 


1359 


1361 


1504 


1658 



The strange baryon masses do not agree perfectly due to our simple model as- 
sumptions. The largest deviation is 41 MeV. The lower bounds on the constituent 
quark mass differences derived from the Feynman-Hellman theorem are fulfilled for 
our fitted quark masses, which are quite close to those of Lichtenbcrg'^^J . The other 
parameters are very similar in value to the usual baryon spectroscopic results^^^l . 

The direct extension of the single baryon Hamiltonian eq.(l) to is neither rea- 
sonable nor successful. First, the two-body confinement interaction will give rise to 
a spurious color van der Waals force. We take this as an indication that the q-q 
confinement interaction between two color singlet hadrons is modified due to the 
nonlinearity of QCD'^^l. The nonperturbative and lattice gauge approaches both give 
rise to a string structure instead of a two-body q-q interaction [-"^^l . The string struc- 
ture and two-body confinement give rise to similar spectroscopic results for simple 
quark system but are not idcntical'^^'^^' . Lattice gauge calculations, after taking the 
(light) qq excitations into account, show that the (heavy) Q — Q confinement inter- 
action is screened. Numerical results can be fitted by the following color screening 
interaction [^^1 

nr)=f-^ + arUi^V (5) 



/ir 

a, = 0.21 ±0.01, 7^^= 400 MeV, = 0.90 ± 0.20 fm. 

Based on these results, we model the Hamiltonian of q^ as an otherwise direct exten- 
sion of eq.(l), but modify the confinement part as 

I — Aj ■ Xjarfj if i,j occur in the same baryon orbit, (qjs^) 

'■^ 1 — Aj • Aj^(l — e'^'^'i) if i,j occur in different baryon orbits. 

Here the exponential e~^^ appearing in eq.(5) has been replaced by a Gaussian e"''^^, 
solely to simplify the numerical calculations. Another reason is this form will auto- 
matically match the quadratic confinement in the short distance {/ir^ <S 1) region. 
Keeping the same form of confinement as that of a single baryon when the interacting 
pair of quarks occur in the same baryon guarantees that when the two baryons are 
separated to large distances, the energy of the q^ system evolves to the two baryon 
internal energy calculated by the Hamiltonian eq.(l). 

B.Svetitskyt^'^] gave a qualitative description of the Q — Q potential: the short- 
range Coulomb potential evolves to a linear part at larger distance, but at still larger 
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distances, the linear part evolves to an exponentially decaying Yukawa potential due 
to light meson exchange. To take this possibility into account, we assume another 
screening confinement potential (as described above). 

_ f ■ ^jO^ffj if hi occur in the same baryon orbit, {'SO) 

\ — Aj • XjarfjC'^^^i if i,j occur in different baryon orbits. 

Although eq.(6A) would be favored as better representing the Yukawa fall-off to a 
constant potential at large distances, that constant is not the zero value expected 
as between colorless hadrons. Thus, eq.(6B) both tests the sensitivity (which turns 
out to be very mild) of our model to the long distance behavior of the potential, and 
provides a model with the expected zero value of the hadronic potential at infinity. 

The screening parameter /i(z/) is determined by fitting our model to the N-N 
scattering phase shifts and they are'^l 

i/ = 0.40fm-2, // = 1.0fm-2. 

Note particularly that the latter is consistent with the lattice result^^^l. 

For the model space, we extend the quark cluster model space by introducing a 
delocalized single quark orbit. In the usual quark cluster model approach, two single 
quark orbits are assumed 

1 \ 4 



0L(r) — \^) ^ (i^f^ centered orbit), 

0ii(r) = ( ^) e (right centered orbit). 



Here s = is the separation of the centers of two cluster. We introduce the 
delocalized quark orbit, 

= (0L(r1+e(s)0«(f))/iV(s), 
^R(r1 = (0jj(f) + e(s)0i(r))/Ar(s), (8) 
N^s) = l + e\s) + 2e{s){<j>L\<l>R). 

The delocalization parameter e(s) is determined variationally for every s by the dy- 
namics of the q^ system (see section III). This orbit is a generalization of the quark 
molecular orbit introduced by Fl.Stancu and L.Wilets^^^l, The six quark space is re- 
stricted to be the space spanned by the following quark cluster bases: 

*aiFi,a.F.(l---6) = ^[^'a.F, (123)^,,^, (456)],, 

V'aiFi(i23) = x{mvs,iM{mM^)M'^)M^) (9) 

^,,^',(456) = x(456)775,/,j,F,(456)V^fl(4)V^,^(5)V^fl(6). 
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Here a = {SIJ) describes the strong interaction conserved quantum numbers of 
strangeness, isospin and spin.Q The cluster WF is almost the same as given in 
eq.(2), but the single cluster Gaussian WF eq.(3b) is replaced by the delocalized 
orbital WF eq.(8), and a flavor symmetry quantum number F is shown explicitly. [ 
refers to isospin and spin coupling by means of the SU2 Clebsch-Gordan coefficients. 
The SU3 color coupling is trivial because only color singlet hadron states are used in 
our calculation. A is the normalized antisymmetry operator 

where p is the two-quark permutation operator. Eq.(9) is termed the physical ba- 
sis by M. Harvey[^°l. Even though only totally symmetric orbital configurations 
are included in eq.(9), q^ orbital excitation configurations are included due to the 
delocalized single quark orbit eq.(8) used in eq.(9). Hidden color channels are not 
included in eq.(9). The reason is that the colorless channels are already complete 
when orbital-spin-isospin excitation configurations are included^^^'^^' . Furthermore, 
color states have not been well constrained in QCD models. It should be clearly kept 
in mind that these physical bases are dependent on the separation s of the two q^ 
cluster centers and the delocalization parameter e(s). 

A check had been done in the N-N scattering dynamical calculation, that if we 
start with the Fl.Stancu and L.Wilets molecular orbits instead of the left and right 
centered orbits (pL and (pR, exactly the same results are obtained'^!. 

III. CALCULATION METHOD 

A dynamical calculation of single channel A^-A^ scattering has been done first 
to fix the screening parameter ^{u) by fitting the model A^-A^ phase shifts to the 
experimental ones'^l 

The q^ states of a given set of quantum numbers a = {SIJ) are expressed as a 
multiple physical channel coupling WF 

V^„(l---6)= C:^F,,a,F,K,Fua,Fr (10) 

aiFi,a2F2 

The channel coupling coefficients p ^^p^ are determined by the diagonalization of 
the Hamiltonian. The maximum number of channels coupled is 16. In the diago- 
nalization, the non-orthogonality property of the physical bases is properly accounted 
for. 

^Orbital angular momentum is assumed to be zero for the lowest states. An angular momentum 
projection which should be done is left for future work. Preliminary estimates indicate this correction 
to the state energy is small. 
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The six quark Hamiltonian matrix elements 

are calculated by the group theory method developed by M. Harvey and by J. Q. 
Chen and ourselves'^"! : 

(1) The physical bases "^aiFi a-zF-z expressed in terms of the symmetry bases 
(group chain classification bases) by the 6—^3x3, SU^„ D SUj„x SU„ and SU3 D 
SUJxUf isoscalar factors calculated by Chen et a/[^^l. 

(2) The six quark Hamiltonian matrix elements in the symmetry bases are reduced 

to a two body matrix element and a four quark overlap (due to the non-orthogonality 
of the delocalized orbit) by the traditional parentage expansion(6 — >■ 4 x 2, SU^n ^ 
SU^xSU^ and SU^ D SUJxUf isoscalar factor[22l). 

(3) The four quark overlap is reduced to a one body overlap by the permutation 
symmetry property of the four quark state^^^l . 

(4) Two body confinement interaction matrix elements are calculated as follows: 
if the interacting quark pair occurs in the same left or right orbit, i.e., {LL\V\LL), 
{RR\V\RR) , {LL\V\RR) and {RR\V\LL) , then the usual quadratic confinement form 
is used. We do so since the two quarks involved arc then always in the same 'baryon', 
and hence in a relative color anti-triplet state. For all the other two body confine- 
ment interaction matrix elements, such as {LR\V\LR), { LL\V\LR), • • •, the color 
screening confinement form is used. Again, we do so since the confining interaction 
will cancel over all the quark pairs identifiable as originating in different color singlet 
hadrons; we effectively remove it in advance from each pairwise interaction leaving an 
interaction which mimics color singlet (mesonic) exchanges. Here {LL\V\RR) means 
(04l)042)|yi2|0ii(l)0fl(2)). 

The eigenenergies obtained in this way are dependent on the separation s and 
the delocalization parameter e(s). We repeat the calculation for each s by varying 
e{s) from 0.1-1.0 with step size 0.1 to get a minimum of the eigenenergy, which thus 
also determines the delocalization parameter e(s). The difference of the minimum 
eigenenergy at separation s and the minimum eigenenergy at infinite separation is 
taken to be the baryon-baryon potential energy Vq{s) (an adiabatic approximation). 
Numerically the asymptotic values of the eigenenergy are indistinguishable from the 
calculation at s = 3fm. It is equal to the threshold sum, i.e., the sum of theoretical 
masses of the corresponding channel baryon pair (single channel case) or of the lightest 
baryon pair (channel coupling case), and the model relative kinetic energy of this 
baryon pair, which is equal to 1/6 of the total kinetic energy of the system due to our 
model WF assumption, eq.(9). (This is one of the checks of our numerical calculation; 
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another check is that all the channel mixing matrix elements are vanishingly small at 
s — 3 fm) . 

A zero point harmonic oscillation energy is added to the minimum potential 

energy 14 (sq) to obtain the binding energy Ba of a state with the quantum number 
a. 

= K(«o) + (11) 

Here /i^ is taken simply to be the reduced mass of the corresponding channel baryon 
pair (single channel case) or the lightest baryon pair within the quantum number a set 
(channel coupling case). In principle, we should do a multichannel coupling dynamical 
calculation; this program is being pursued only for the few most promising dibaryon 
candidates. 

Finally the experimental channel baryon pair mass (single channel case) or the 
lightest baryon pair mass (Mi+M2)a (channel coupling case) is added to Bo. to obtain 
the lowest mass of each quantum number a set, 



M,(g6) = (Mi + M2)« + S,. (12) 

The mass Ma{q^) is compared not only to the two body decay threshold, the experi- 
mental hghtest baryon pair mass (Mi + M2)a, but also to the possible multi-particle 
final states allowed by strong interaction to determine if there is a strong interaction 
quasi-stable dibaryon state. 

In order to show the effects of channel coupling, flavor-symmetry breaking 

and the different forms of color screening, the following eight sets of calculations 
have been done: 

(1) single channel, flavor symmetry (scs), 

(2) multi-channel, flavor symmetry (ccs), 

(3) single channel, flavor symmetry breaking (scb), 

(4) multi-channel, flavor symmetry breaking (ccb), 

where (l)-(4) have been calculated using the color screening form (6A) and (5)- (8) 
replace this with the color screening form (6B). 

To indicate the level of uncertainty due to the choice of the color screening pa- 
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rameter value, we have calculated the results corresponding to n = 1.6/m~^ and 
u = 0.46, 0.60/m~^ in addition to the best fit values, = 1.0/m~^ and u = 0.40/m~^. 
(The reasoning behind the parameter values chosen for these additional cases is de- 
scribed in the Appendix.) Only dibaryon candidates based on both the channel 
coupling and symmetry breaking results should be considered serious possibilities for 
experimental searches. 

A computer program package which incorporates all the needed group theory 
results has been written to automate the numerical calculation. It can be used for 
other model calculations by simply replacing the one- and two-body matrix elements. 
In particular, it is may also be used for a relativistic quark cluster model calculation. 
As a cross-check on the program, the S = 0, (IJ) = (01), (10) and (03) channels have 
been done both by the method described above and by direct diagonalization with 
the physical bases. 

IV. RESULTS 

All possible sets oi a — (SIJ) within the u, d, and s three flavor world have been 
calculated. Only a few states are strong interaction quasi-stable or narrow resonances, 
and they are hsted in Table II. 

The strong interaction unstable states have been omitted to simplify the presen- 
tation. However some general features are listed here. 

(1) The two color screening forms given quite similar results: The form (6A) gives 
slightly higher six quark masses than the form (6B) but the largest difference is only 
about 10 MeV for the {SIJ) = (003) case. This is consistent with the findings for 
N-N scattering where the forms (6A) and (6B) give similar N-N phase shifts while 
form (6B) yields a httle stronger attraction in the ^5"! channel ^1. 

(2) There are two extreme kinds of dibaryon candidates. One kind is a loosely 
bound ('molecular') two baryon state. The binding energy is small (usually around 
zero), J < 1; the delocalization e(so) is also small (usually ~ 0.2). (The deuteron is 
a typical example while the H particle is an exception.) Their masses are close to 
the lowest two body decay threshold and they might be stable with respect to the 
strong interactions, but their stability is very sensitive to model details. The other 
kind is a tightly bound state. The binding energy is large (~ 100 MeV), J > 2; 
the delocalization is also large (~ 0.8). Their masses are larger than the lowest two 
body threshold and they are therefore unstable. However their masses are lower 
than the favorable three or four body decay threshold. Their two or three body 
decay is hindered due to large angular momentum and so they might nonetheless 
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appear as narrow resonances. A typical example is the di-A(003) state. This quasi- 
stability property is not sensitive to the model details, unless radically different model 
assumptions are made. 

(3) Flavor symmetry breaking effects are channel dependent as reported by Malt- 
mant^^], decreasing the binding by an amount ranging from nil to 70MeV. For some 
multi-channel coupling cases, the flavor symmetry breaking effect changes which is 
the lowest channel. In those cases, there is large apparent flavor symmetry breaking 
effect {ABa ~ 120 — ISOMeV). The mass of the state is increased correspondingly 
and this produces a large difference in the stability of the state with respect to the 
multi-particle decay channels, due to the flavor symmetry breaking. 

(4) The channel couphng effect is small for most cases, even after taking into 
account quark delocalization in the extended u, d, s world. For the H particle case, 
channel crossing occurs: The dominant channel for the minimum potential (SS) and 
the asymptotic channel (AA) are different. This channel coupling decreases the q^ 
mass by about 100 MeV and makes the H particle stable. A similar channel coupling 
effect occurs in the (S'/J) = (-202) state. 

Wc emphasize that all the QDCSM parameters are fixed by the ground state 
baryon properties and the N-N scattering phase shifts. Therefore the dibaryon can- 
didates noted here are a relatively robust theoretical prediction. We wish to emphasize 
a few additional points. 

(1) In the S = 0, (/J) = (01) channel (the deuteron channel), the model predicts 
that there is a state with M(OOl) = 1880(±14) MeV. It is a di-nucleon state because 
the delocalization e = 0.2 is small and 4 MeV away from the deuteron energy. We 
take this as a measure of the predictive power of QDCSM for the dibaryon state. Due 
to this uncertainty, wc have included those states which are close to the lowest two 
body decay threshold in table II. 

(2) In the -S" = 0, (/J) = (10) channel, the model predicts that there is a di- 
nucleon (e(so) = 0.1) resonance state at 1889(±4) MeV. It is about 10 MeV away 
from a possible zero binding di-nuclcon resonance. This is another example that 
shows the predictive power of QDCSM and appears to limit the model uncertainty. 

(3) In the nonstrange sector (S — 0), the model predicts an a = (003) state with 
mass M(003) = 2110(±36) MeV. It has the largest binding (320-390 MeV) within the 
u, d, s three flavor world. Although its mass is above NNn threshold, the transition 
to NNtt and NN is hindered by the large angular momentum, so that it is still 
possible for this to be a narrow resonance. The large delocahzation, e = 1.0, means 
this is a true six quark state. All these results are consistent with our earlier simple 
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relativistic model result'^^^. Although the skyrmion model calculation of Walet does 
not obtain a binding as large as QDCSM result'^^^, we note that this model doesn't 
obtain sufficient attraction in the N-N channel, either. We know of no reason for the 
predictive power of the QDCSM shown in the N-N channel to be totally lost in the A- 
A channel since both channel results arc stable under the same reasonable variations 
in the values of the model parameters. Therefore, we continue to recommend this 
state highly as a good candidate for discovery of a dibaryon resonance. 

There is no other interesting channel in the u, d two flavor world as found earlier 
by Maltmanl'^l . 

(4) Jaffe's H particle remains as the unique strong interaction stable dibaryon 
in QDCSM. M{H) = 2199(±24) MeV is 32 MeV lower than the AA threshold. The 
delocalization e(so) = 1.0 is also large and the adiabatic channel coupling WF is quite 
close to Jaffe's pure symmetric flavor singlet basis. However, due to the sensitivity to 
details of the model, it is not possible to claim that it is indeed a strong interaction 
stable state. 

(5) Another interesting state is the a = (-3|2) statct^^l Af(-3|2) = 2529(±25) 
MeV is 45 MeV lower than the favorable (AStt) three body decay threshold. Also 
due to large angular momentum, its decay into AS should be inhibited and so it too 
might show up as a dibaryon resonance. Another good point about this state is that 
all the other states with the same quantum number set are about 100 MeV higher 
than it (not listed in table II). This might make it a cleaner resonance to observe. 

(6) The states M(-l|0) = 2133(±7) MeV,M(-220) = 2390(±9) MeV,M(-3|l) 
= 2512(±15) MeV, M(-3il) = 2469(±58) MeV, M(-400) = 2637(±10) MeV,M 
(—600) = 3359(±18) MeV, all have their masses close to the corresponding thresh- 
olds: NE, EE, ES, AS, SS and nn. Another group of states M(-li3) = 2318(±33) 
MeV, Af(-213) = 2530(±37) MeV, M(-202) = 2345(±26) MeV, M(-3|3) = 2728 
(±37) MeV, have very large binding and their masses are all less than the favorable 
multi-body channel. These states bear further study. The high spin dibaryon res- 
onances seem to be especially worth experimental searches, in addition to the spin 
zero H particle. 

V. CONCLUSION 

Since Jaffe's first prediction of H particle^^l, there have been many efforts both 
theoretically and experimentally to search for dibaryons. Whether all these QCD in- 
spired models miss some physics when they are extended from the single hadron to the 
multi-hadron case, so that their predictions are not reliable, remains a question. The 
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QDCSM is a more realistic model by taking into account the possible difference of the 
q — q interaction inside a single baryon and between two color singlet baryons, by al- 
lowing each system to choose its favorable configuration in a larger Hilbert space, and 
by having the model constrained not only by qualitatively fitting hadron spectroscopy 
but also N-N scattering. This model approach has some moderate success to sup- 
port it: It predicts two di-nucleon states not too far from the experimentally known 
deuteron and quasi-deuteron states. If we take this as a measure of the predictive 
accuracy of the QDCSM, then there are few promising dibaryon resonance candidates 
within the u, d and s three flavor world, as listed in table II. Due to the simplicity 
of the model assumptions, the quantitative predictions of the dibaryon masses are 
uncertain (~ lOMeV for nonstrange states and even larger uncertainty for strange 
states). To get more reliable estimates, especially to be able to determine whether or 
not the candidate states are strong interaction stable, requires improvement of many 
aspects of the QDCSM: 

(1) The QDCSM does not fit the baryon octet and decuplet perfectly, the largest 
deviation being 41 MeV for the S. Although the adiabatic potential is obtained 
through a subtraction procedure which suggests cancellation of errors is possible, 
there is no guarantee that the uncertainty of the strange baryon mass cancels very 
accurately. 

(2) An adiabatic approximation has been used in this calculation which should be 
replaced by a dynamical channel coupling calculation. 

(3) This calculation is nonrelativistic; a relativistic calculation is underway to 
estimate relativistic corrections. The preliminary result is that the relativistic and 
nonrelativistic versions give very similar mass values, especially for the nonstrange 
states. 

(4) Only A^-A^ scattering has been used to constrain the QDCSM. Although data 
is sparser, A-p and S-p scattering should be used as well. We have begun such an 
analysis. 

(5) The effects of qq excitations or quark-meson couplings, which have not been 
included, may well be important, especially for those states which have a mass close 
to the lowest two-body threshold. 

(6) It would be interesting to include c and b quarks with a view towards making 
contact with Heavy Quark Effective Theory. However, the large quark WF difference 
between u, d and b, c would need to be treated first. This is unlike the s quark case 
where the single quark WF distortion is not large. 
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Wc believe the QDCSM results support the value of investing additional effort, 
both theoretical and experimental, using more sophisticated approaches, in order to 
concentrate on a few promising dibaryon candidates. 
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APPENDIX 



In the dynamical N-N scattering calculation PI , we first calculate the interaction 
kernel 

KCs.s)-^m!m^ (An 

< ■^{s)\-^{s) > 

where H is the six quark Hamiltonian, "^{S) is the N-N channel WF (9). Then we 
do a partial wave decomposition 

ki{S, S') = J dnPi{Qss')K{S, S') {A2) 

where Qgs' means the angle between S and S'. We then assume the diagonal matrix 
elements ki{s, s) as the effective interaction of the 1-partial wave between two nucleons, 

Vi{s) = ki{s,s) {A3) 

This effective interaction is dependent on the delocalization parameter e(s). We vary 
the value e(s) to get a minimum for each separation s. The e(s) so determined is I 
dependent, i.e., we have ei{s). Next we substitute the values ei{s) and e;(s') back into 
eq. (A2) to get the final ki{s, s') for the N-N scattering calculation. Finally, wc adjust 
the screening parameter or z/ of cq. (6 A) and (6b) to get the best fit to the ^5*0 and 
^Si N-N phase shifts. This determines the best values: = 1.0/m~^, u = 0.4/m~^, 
because it is a more complete and consistent calculation. However it takes much more 
computer time than the next calculation we describe. 

To minimize computer time for systematic dibaryon search, we tried another ap- 
proximation, assuming K{s, s) as the effective interaction between two nucleons 

V{s) = K(s, s) (^4) 

Then, we varied e(s) to minimize V{s) and so determine the e(s) as well. This e(s) is 
partial wave independent. Substituting the e(s) and e(s') values so determined back 

to eq. (Al), we obtain the final K{s,s') and then do a partial wave decomposition 
to calculate the phase shifts. Adjusting the screening parameters again to obtain the 
best fit, the ^5*0 and ^Si channels, we obtain a second set of values: /i — 1.6/m~^, 
u — 0.6/m~^. The approximation is not as good as the first one, but also gives a 
quafitatively good fit to the N-N phase shifts. 

In the dibaryon calculation, we use the second variation method (variation before 
partial wave decomposition) to obtain the effective interaction between two baryons. 
The results are shown in Tables Ila and lie. As a check on the range of variation, we 
also used the intermediate value u — 0.46/m~^, and those results are shown in Table 
lib. 
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Table Ila. = 1.0/m-^ u = 0.4/m-^ 
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Table lib, v = 0.46/m-^ 
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Table lie, jj, = 1.6 fm ^, = 0.6/m ^ 
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2 


0.2 


1.4 




-200 


scs 


2097 


-244 


-161 


1.0 


0.8 


2097 


-244 


-161 


0.9 


0.8 


2231 (AA) 




CCS 


1974 


-364 


-257 


1.0 


0.7 


2012 


-365 


-220 


1.0 


0.6 






scb 


2282 


-240 


-104 


1.0 


0.6 


2281 


-241 


-105 


1.0 


0.6 






ccb 


2176 


-200 


-55 


1.0 


0.6 


2175 


-202 


-56 


1.0 


0.6 




-202 


scs 


2242 


-272 


-231 


0.8 


1.1 


2249 


-265 


-223 


0.7 


1.1 


2472(NS*) 




CCS 


2150 


-373 


-323 


1.0 


1.0 


2255 


-368 


-318 


1.0 


1.0 


2397(NS7r) 




scb 


2413 


-221 


-164 


0.9 


0.9 


2429 


-219 


-148 


1.0 


0.8 






ccb 


2319 


-231 


-153 


1.0 


0.8 


2320 


-230 


-152 


1.0 


0.8 




-213 


scs 


2372 


-423 


-393 


1.0 


1.2 


2382 


-413 


-383 


1.0 


1.2 


2765(AS*) 




CCS 


2372 


-423 


-393 


1.0 


1.2 


2382 


-413 


-383 


1.0 


1.2 


2690(AS7r) 




scb 


2496 


-308 


-273 


1.0 


1.1 


2509 


-302 


-260 


1.0 


1.0 


2511(AA7r7r) 




ccb 


2493 


-308 


-272 


1.0 


1.1 


2506 


-302 


-259 


1.0 


1.0 




-220 


scs 


2391 


-20 


5 


0.2 


1.4 


2388 


-18 


2 


0.1 


1.6 


2386(SS) 




CCS 


2390 


-21 


4 


0.2 


1.4 


2388 


-18 


2 


0.1 


1.6 






scb 


2382 


-33 


-4 


0.2 


1.3 


2386 


-29 





0.2 


1.3 






ccb 


2381 


-34 


-5 


0.2 


1.3 


2385 


-30 


-1 


0.2 


1.3 




-3|3 


scs 


2510 


-423 


-394 


1.0 


1.2 


2520 


-413 


-384 


1.0 


1.2 


2904(A17) 




CCS 


2510 


-423 


-394 


1.0 


1.2 


2520 


-413 


-384 


1.0 


1.2 


2788(AS*7r) 




scb 


2704 


-254 


-214 


1.0 


1.0 


2717 


-250 


-201 


1.0 


0.9 


2714(AS7r7r) 




ccb 


2691 


-255 


-213 


1.0 


1.0 


2705 


-250 


-200 


1.0 


0.9 




-3fl 


scs 


2493 


-46 


-18 


0.2 


1.3 


2497 


-42 


14 


0.2 


1.3 


2511(SS) 




CCS 


2491 


-48 


-20 


0.3 


1.3 


2497 


-42 


14 


0.2 


1.3 






scb 


2498 


-60 


13 


0.4 


1.0 


2501 


-57 


10 


0.3 


1.0 






ccb 


2497 


-61 


14 


0.4 


1.0 


2500 


-58 


11 


0.4 


1.0 




-3i2 

Z 


scs 


2339 


-321 


-272 


1.0 


1.0 


2344 


-316 


-268 


1.0 


1.0 


2611(Nf7) 




CCS 


2287 


-373 


-324 


1.0 


1.0 


2292 


-368 


-319 


1.0 


1.0 


2574(AS7r) 




scb 


2510 


-313 


-193 


1.0 


0.6 


2508 


-304 


-194 


1.0 


0.6 






ccb 


2506 


-240 


-106 


1.0 


0.6 


2504 


-242 


-107 


1.0 


0.6 




-3il 


scs 


2339 


-321 


-272 


1.0 


1.0 


2344 


-316 


-268 


1.0 


1.0 


2434(A5) 




CCS 


2337 


-156 


-97 


1.0 


0.9 


2340 


-154 


-94 


1.0 


0.9 






scb 


2474 


-167 


-38 


1.0 


0.6 


2437 


-131 


3 


1.0 


0.6 






ccb 


2422 


-146 


11 


1.0 


0.6 


2421 


-147 


13 


1.0 


0.6 




-400 


scs 


2617 


-41 


-19 


0.2 


1.4 


2622 


-37 


-15 


0.2 


1.4 


2636(SS) 




CCS 


2617 


-41 


-19 


0.2 


1.4 


2622 


-37 


-15 


0.2 


1.4 






scb 


2628 


-45 


-8 


0.3 


1.1 


2630 


-43 


-6 


0.2 


1.1 






ccb 


2627 


-46 


-9 


0.3 


1.1 


2630 


-43 


-6 


0.2 


1.1 




-600 


scs 


3251 


-112 


-94 


0.3 


1.4 


3260 


-102 


-84 


0.3 


1.4 


3345(Of2) 




scb 


3376 


-65 


32 


1.0 


0.6 


3375 


-67 


30 


1.0 


0.6 
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